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» I 1 We present analytical and numerical studies of the Fourier transform (FT) of the 

' gravitational wave (GW) signal from a pulsar, taking into account the rotation of 

the Earth for a one day observation period. 

1 INTRODUCTION 

The direct detection of gravitational radiation (GR) from astrophysical sources is 
one of the most important outstanding problems in experimental gravitation today. 
The constructiorL.of large laser interferometric gravitational wave detectors like the 
LIGOly, VIRGO H, LISA, TAMA 300, GITO 6000 and AIGO is opening a new win- 
dow for the study of a vast and rich variety of nonlinear curvature phenomena. The 
network of gravitational wave detectors can confirm that GW exist and by mon- 
itoring gravitational wave forms give important information on their amplitudes, 
frequencies and other important physical parameters. 

A prototype of continuous astrophysical sources is a pulsar. A variety of in- 
stabilities cause deformations from spherical symmetry giving rise to GWs. The 
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amplitudes of GR from these pulsars are probably very weak(< 10 -26 — 10~ 28 , for 
galactic pulsars). The GR signal will be buried deep within the noise of the detec- 
tor system. The detection of a GR signal warrants the urgent need of careful data 
analysis with development of analytical methods and problem oriented algorithms. 

In sections 2, 3 and 4 we outline the approach leading to the FT of the GW 
signal. The frequency modulation (FM), Doppler shift due to rotation and orbital 
motion of the Earth in the Solar System Barycentre (SSB) frame, its effect on the 
total phase of the received GW signal and the Fourier transform (FT) of the GW 
signal have previously been describedu. 

Sections 5 and 6 present discussion and conclusions. 

2 Methodology 

Typical values of the gravitational wave amplitude h for the Crab and Vela pulsars 
are ~ 10 -25 and ~ 10~ 24 respectively. This amplitude is several orders of magni- 
tude below LIGO's expected sensitivity of ~ 10 -23 ). Since the LIGO would make 
continuous observations over a time scale of a few months or more, a significant 
enhancement to the signal-to-noise ratio (SNR) is expected by integrating the data 
over a long time interval. 

The total response of the detector is a function of the source position, the 
detector orientation, the orientation of the spin axis of the Earth and the orientation 
of the orbital plane. Since the pulsar signal is weak, long integration times ~ 10 7 
sees will be needed to extract the signal from the noise. Since the detector moves 
along with the Earth in this time, the frequency of the wave emitted by the source 
is Doppler shifted. Also since the detector has an anisotropic response, the signal 
recorded by the detector is both frequency and amplitude modulated. We discuss 
now the important role of frequency modulation in the context of signal detection. 

3 Study of Frequency Modulated Pulsar Signal 

Frequency modulation arises due to translatory motion of the detector acquired 
from the motion of the Earth. We consider only two motions of the Earth: its 
rotation about the spin axis and the orbital motion about the Sun, so the response 
is doubly frequency modulated with one period corresponding to a day and the 
other period corresponding to a year. The FM smears out a monochromatic signal 
into a small bandwidth around the signal frequency of the monochromatic waves. It 
also redistributes the power in a small bandwidth. The study of FM due to rotation 
of the Earth about its spin axis for a one day observation period shows that the 
Doppler spread in the angular bandwidth for 1kHz signal will be 0.029Hz. The 
Doppler spread in the angular bandwidth due to orbital motion for an observation 
period of one day will be 1.7 x 10 _3 Hz 0. 

Since any observation is likely to last longer than a day, it is important to 
incorporate this effect in the data analysis algorithms. 

In order to study frequency modulation of a monochromatic plane wave, one 
must calculate the Doppler shift due to rotational and orbital motion of the Earth 
in the SSB frame. For this, we need to know relative velocity between the source 
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and detector. The Euler angles (6, <j>) give the direction of the incoming wave in 
the SSB frame. We characterize the motion of the Earth (and detector) simply: 
(a) We assume the orbit of the Earth to be circular, (b) We neglect the effect of 
the Moon on the motion of the Earth. 

The phase <f>(t) of the Doppler shifted received signal for a single direction sky 
search (#,</>) is given by, 



4>(t) = 2tt / f rec (f) dt' 
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where 0' = Lo orh t - cj>, /3' = /3 + Wroti, O = ^or&io - 4>, P' = (3 + uj rot t 0l [3 
is the initial azimuthal angle of the detector at the observation time to- where 
A is distance from the centre of the SSB frame to the centre of the Earth, R 
is the radius of the Earth, and n the unit vector in the direction of source, 
n = (sin#cos</>, sin 9 sin <f>, cos#). Here we have assumed that at time t = 0, 
the longitudinal angle (3 = 0. 

It can be seen from eq (||) that the Doppler corrections to the phase of received 
pulsar signal depends on the direction of the source in the sky. Orbital motion is 
omitted in this preliminary analysis but will be presented in later work. 



4 Fourier Transform Analysis of the FM signal due to the 
Rotational Motion of the Earth 

We analyze the Fourier transform (FT) of the frequency modulated signal and 
study the extent to which the peak of the FT is smudged and to how much the FT 
spreads in the frequency space. This type of study would be useful from the point 
of view pi data analysis and for applying such schemes as stepping around the sky 
method □ which relies on the FT. 

We let x = 2 ^fc> R and t = 0. Here x plays the role of modulation index similar 
to K in the theory of signal modulation. The modulation index depends on the 
frequency of the pulsar signal. If we consider only the frequency modulated output 
of the signal, the output of amplitude unity is given as follows, 

h{t) = cos (</>(£)) . (4) 

We now consider the h (t) to be given on a finite time interval [0, T] which would 
be assumed to be the observation period. In our analysis we have assumed T to be 
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one day. The Fourier transform of the signal h (t) is given by, 



h(f) = f h{t)e- l2nft dt. (5) 



(i 



It is convenient to use a time coordinate £ = uj ro tt/2 which for a period of a day is 
of the order of unity, i.e. when T = lday = 86400 sees, the £ T = w rot T = ir. 

An exact closed form expression for the Fourier transform of the frequency 
modulated GW signal is obtained by the- analytical approach. The plane wave 
expansion in spherical harmonics is used, u 

Then the GW signal becomes 



^2 J Si m (fo,a,t)dt 



An[Y / ^A0A)e 27rlfot Ji(ksm(cj r t/2))Y lm (a, -L-)]dt (6) 
P =0 1 

where K is the wave vector with spherical coordinates (K, 6, <f>) with 

\K\ = 4tt/ 5R sin a(wa.(u r t/2)) = k sin(w r t/2) (7) 

and h = (n, 0, (j>) (n = 1). 
and the FT is given by 

Simif, fo, a) = FT(S lm (f , a, t)) (8) 
Further simpification gives: 

lm I m 

where 
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and 
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With convenient co-ordinate transformations, we obtain 

Si m (jj,wo,a) = Ai m I(l,B, k) (12) 

where 



I(l,B,k) = exp[iBQ ' — d£_ 

Jo y/ksmt, 

1(1, B,k) = 2i B [ 
Jo 



Also, 

^ cos(.B cos -1 x) Ji + i/ 2 {kx) 



kx\f\ 



(13) 



dx (14) 



astro2: submitted to World Scientific on February 1, 2008 



(where B = 2{^f-) +m ) 
which equals 



2tt 



/ " k l+1 / 2 T(l + 1) iF 3 (l + l;l + 3/2, (Z + 2 + B)/2, (Z + 2 - B)/2; -fc 2 /16) 



^2 2 2 '+! r(Z + 3/2) I\(Z + 2 + S)/2)r((i + 2 - fl)/2) 

(15) 

The 1-F3 is the generalised hypergeometric function with one parameter (here 
I + 1) in the numerator and three parameters in the denominator and — k 2 /16 as 
its variable. A wealth of information is hidden in these parameters and 1-F3 is 
uniformly convergent. 

The numerical integration is done via the Clcnshaw-Curtis quadrature algorithm 
and gives accurate and reliable results. Both the analytical and numerical results 
agree for k < 1 and arbitrary /; a result of possible relevance for LISA which 
observes the low frequency band 0.1 to 0.00001 Hz. For simultaneously large k and 
I, the factor feC+V 2 ) presents a challenge. 

We also evaluate 

I c tt fc 2 ' r(l/2)r(Z + l/2) 1 F 2 (Z + l/2;2Z + 2,Z + 3/2;-fc 2 ) 

1 (mWI 4 2 2 <+! r(Z + l)r(Z + 3/2)r(Z + 3/2) [ ' 

In this case we find perfect agreement with numerical results for k < 86 and ar- 
bitrary Z(< 300). The above approach using the plane wave representation appears 
to be equally useful when orbital motion is also included in the analysis and will 
be discussed in future work involving orbital corrections. 



5 Discussion 



In general, the FT is a complex valued function of the parameters k, I and B. 
Where appropriate, the absolute, real and the imaginary parts are taken. In Figure 
1, the absolute value of the integral is plotted as a function of k for B = 30 and 
I = 10. The maximum occurs around k = 40, which translates to approximately 
150 Hz. The intensity drops off fast and is negligible around k = 300. Figure 2 
again shows the absolute value of I vs k for B = 10 but for I = 270 >> 1. This 
plot shows remarkable evidence for the Bessel function index theorem as it shows 
I(k) = for k < 250, rising to a maximum around k — 285 and falling off fast 
thereafter. For both large k and Z, B must be small to control the oscillations of 
the integrand to give a non-negligible value. Figure 3 shows absolute I vs B for 
I = k = 10. The integral peaks at B = and drops off fast around B = 6.25 
and rises again to a tiny value before dropping off to zero B > 10. Figure 4 shows 
absolute I vs B for I = k = 200(>> 1). The integral peaks at B = then drops 
off rapidly at B > 20. This again illustrates that B has to be small to control the 
oscillations. Figure 5 shows absolute / vs B for Z = 10 and k = 270(>> 1). Here 
the many bands of B are in full display until B < 280. Figure 6 shows the absolute 
value \Si m (k)\ 2 vs Z and k for the indicated range. 
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Figure 1. |J(Z, fc, B)| vs k for B = 30, 1 = 10 



6 Conclusion 

We have studied the frequency modulation of the pulsar signal by two different 
methods. A closed form of the Fourier transform of the frequency modulated GW 
pulsar signal due to rotational motion of the Earth about its spin axis has been 
obtained. The work with the inclusion of orbital corrections results in a double 
series with only a finite number of surviving terms.. Further details of this analysis 
warrant study and will be presented in later workfi Deeper results would be useful 
for schemes like the stepping around the sky methodB and for differential geometric 
methods that allow the setup of search templates in the relevant parameter spaceu. 
This analysis should also be applicable to the Doppler modulation of the GW signal 
caused by LISA's orbit around the sun. 

References 

1. A. Abramovici et al, Science, 256, 325 (1992). 

2. K. Jotania, S. Valluri, and S.V. Dhurandhar, Astron & Astrophys 306, 317 
(1996). 

3. B.F. Schutz, in The Detection of Gravitational Waves, ed. D.G. Blair (Cam- 
bridge University Press, 1991). 

4. B.S. Sathyaprakash, Filtering gravitational waves from Supermassive Black 
Hole Binaries, LISA Conference, CALTECH(1998). 

5. B.H. Bransden & C. J. Joachain, Intoduction to Quantum Mechanics, John 



astro2: submitted to World Scientific on February 1, 2008 6 



Figure 2. \I(l, k, B)\ vs k for B = 10,1 = 270, width=10cm, height=9cm 
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